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Abstract 

Lyapunov exponents of a dynamical system are a useful tool to 
gauge the stability and complexity of the system. This paper of- 
fers a definition of Lyapunov exponents for a sequence of free linear 
operators. The definition is based on the concept of the extended 
Fuglede-Kadison determinant. We establish the existence of Lya- 
punov exponents, derive formulas for their calculation, and show 
that Lyapunov exponents of free variables are additive with respect 
to operator product. We illustrate these results using an exam- 
ple of free operators whose singular values are distributed by the 
Marchenko-Pastur law, and relate this example to C. M. Newman's 
"triangle" law for the distribution of Lyapunov exponents of large 
random matrices with independent Gaussian entries. As an interest- 
ing by-product of our results, we derive a relation between the ex- 
tended Fuglede-Kadison determinant and Voiculescu's S-transform. 

1 Introduction 

Suppose that at each moment of time, £j, a system is described by a state 
function ip (tj) and evolves according to the law ip (tj+i) = X^ip (ij) , where Aj 
is a sequence of linear operators. One can ask how small changes in the initial 
position of the system are reflected in its long-term behavior. If operators Xj 
do not depend on time, Xj = X, then the long-term behavior depends to a large 
extent on the spectrum of the operator X. If operators Xi do depend on time but 
can be modelled as a stationary stochastic process, then the long-term behavior 
of the system depends to a large extent on so-called Lyapunov exponents of the 
process A,. 

The largest Lyapunov exponent of a sequence of random matrices was in- 
vestigated in a pioneering paper Q by Furstenberg and Kesten. This study 
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was followed in [13| by Oseledec, who researched other Lyapunov exponents 
and finer aspects of the asymptotic behavior of matrix products. These inves- 
tigations were greatly expanded and clarified by many other researchers. In 



particular, in LJ] Ruelle developed a theory of Lyapunov exponents for random 
compact linear operators acting on a Hilbert space. Lyapunov exponents for 
random N x N matrices when N — > oo were studied in 0], 10], and 0- 

The goal of this paper is to investigate how the concept of Lyapunov expo- 
nents can be extended to the case of free linear operators. It was noted recently 



by Voiculescu ([17|) that the theory of free operators can be a natural asymp- 
totic approximation for the theory of large random matrices. Moreover, it was 
noted that certain difficult calculations from the theory of large random ma- 
trices become significantly simpler if similar calculations are performed using 
free operators. For this reason it is interesting to study whether the concept 
of Lyapunov exponents is extendable to free operators, and what methods for 
calculation of Lyapunov exponents are available in this setting. 

Free operators are not random in the traditional sense so the usual definition 
of Lyapunov exponents cannot be applied directly. Our definition of Lyapunov 
exponents is based on the observation that in the case of random matrices, 
the sum of logarithms of the k largest Lyapunov exponents equals the rate at 
which a random fc-dimensional volume element grows asymptotically when we 
consecutively apply operators Xi. 

In the case of free operators we employ the same idea. However, in this 
case we have to clarify how to measure the change in the " ^-dimensional volume 
element" after we apply operators Xi. It turns out that we can measure this 
change by a suitable extension of the Fuglede-Kadison determinant. Given this 
extension, the definition proceeds as follows: Take a subspace of the Hilbert 
space, such that the corresponding projection is free from all and has the 
dimension t relative to the given trace. Next, act on this subspace by the 
sequence of free operators X^. Apply the determinant to measure how the 
"volume element" in this subspace changes under these linear transformations. 
Use the asymptotic growth in the determinant to define the Lyapunov exponent 
corresponding to the dimension t. 

It turns out that the growth of the i-dimensional volume element is expo- 
nential with a rate which is a function of the dimension t. We call this rate the 
integrated Lyapunov exponent. It is an analogue of the sum of the k largest Lya- 
punov exponents in the finite-dimensional case. The derivative of this function 
is called the marginal Lyapunov exponent. Its value at a point t is an analogue 
of the fc-th largest Lyapunov exponent. 

Next, we relate the marginal Lyapunov exponent fx (t) to the Voiculescu 
^-transform of the random variable X* Xi. The relationship is very simple: 

fx {t) = -{l/2)log[S x . x {-€)]. (1) 
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Using this formula, we prove that the marginal Lyapunov exponent is de- 
creasing in t, and derive an expression for the largest Lyapunov exponent. For- 
mula ([T]) also allows us to prove the additivity of the marginal Lyapunov ex- 
ponent with respect to operator product: If X and Y are free, then fxy (i) = 
fx (t) + h (t) . 

As an example of application of formula ([1]), we calculate Lyapunov expo- 
nents for variables X; such that X*Xi are distributed as Free Poisson variables 
with parameter A. The case A = 1 corresponds to the random matrix case con- 
sidered by C. M. Newman in ]9l, and the results of this paper are in agreement 
with Newman's "triangle" law. In addition, our results regarding the largest 
Lyapunov exponent agree with the results regarding the norm of products of 
large random matrices in Q . Finally, our formula for computation of Lyapunov 
exponents seems to be easier to apply than the non-linear integral transforma- 
tion developed in 0]. 

An interesting by-product of our results is a relation between the extended 
Fuglede-Kadison determinant and the Voiculescu S'-transform, which allows ex- 
pressing each of them in terms of the other. In particular, if Y is a positive 
operator and if {Pt} is a family of projections which are free of Y and such that 
E(P t ) = t, then 

logSy (-*) = -2j t [logdet {^YPt)\ • (2) 
and if X is bounded and invertible, then 

i r 1 

logdet (X) = --/ log S x *x{-t)dt. (3) 
* Jo 

Calculations related to © and (O were performed by Haagerup and Larsen 
in Q in their investigation of the Brown measure of i?-diagonal operators. The 
Brown measure of an operator X is closely related to the determinant of X — 
zl, and Haagerup and Larsen computed the Brown measure of an i?-diagonal 
operator X in terms of the S-transform of X*X. However, it appears that 
formulas @ and ^ have not been stated explicitly in 



In addition, Sniady and Speicher showed in 15( that an i?-diagonal X can be 
represented in the triangular form and that the spectra of the diagonal elements 
in this representation satisfy certain inequalities in terms of the S'-transform of 
X*X. Sniady and Speicher used their result to give a different proof of the 
results in [5|]. It is likely that Sniady and Speicher's method can also be used 
for a different proof of formulas ^ and ([3]). 

The rest of the paper is organized as follows: Section [2] describes the exten- 
sion of the Fuglede-Kadison determinant that we use in this paper. Section [3] 
defines the Lyapunov exponents of free operators, proves an existence theorem, 
and derives a formula for the calculation of Lyapunov exponents. Section [?] 
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computes the Lyapunov exponents for a particular example. Section [5] connects 
the marginal Lyapunov exponents and the S'-transform, proves additivity and 
monotonicity of the marginal Lyapunov exponent, and derives a formula for 
the largest Lyapunov exponent. In addition, it derives a relation between the 
determinant and the S'-transform. And Section [5] concludes. 

2 A modification of the Fuglede-Kadison deter- 
minant 

Let A be a finite von Neumann algebra and E be a trace in this algebra. Recall 
that if X is an element of A that has a bounded inverse, then the Fuglede- 
Kadison determinant (Q) is defined by the following formula: 

det (X) = exp -E log (X*X) . (4) 

The most important property of the Fuglede-Kadison determinant is its multi- 
plicativity: 

det (XY) = det (X) det (Y) . (5) 

This determinant cannot be extended (non-trivially) to operators with non-zero 
kernel if we require that property ((5]) holds for all X and Y. 

However, if we do not insist on this property, then we can define an extended 
determinant as follows: Let log +A (£) =: logt if t > A and =: if t < A. Note 
that i?log +A (X*X) is a (weakly) decreasing function of A on the interval (0, 1) , 
and therefore it converges to a limit (possibly infinite) as A — * 0. 

Definition 1 

det {X) = exp i lim£:iog +A {X* X) . 

This extension of the Fuglede-Kadison determinant coincides with the ex- 
tension introduced in Section 3.2 of Q . 

Example 1 Zero operator 

From Definition [U if X = 0, then det X = 1. 

Example 2 Finite dimensional operator 

Consider the algebra of n-by-n matrices M n (C) with the trace given as the 
normalization of the usual matrix trace: E (X) = n _1 Tr (X) . Then the original 
Fuglede-Kadison determinant is defined for all full-rank matrices and equals the 
product of the singular values of the operator in the power of 1/n. It is easy 
to see that this equals the absolute value of the usual matrix determinant in 
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the power of l/n. The extended Fuglede-Kadison determinant is defined for 
all matrices, including the matrices of rank k < n, and equals the product of 
non-zero singular values in the power of l/n. 

We can write the definition of the determinant in a slightly different form. 
Recall that for a self-adjoint operator X £ A we can define its spectral probability 
measure as follows: 



where S is an arbitrary Borel- measurable set and lg is its indicator function. 
Then, the determinant of operator X can be written as 



For all invertible X the extended determinant defines the same object as the 
usual Fuglede-Kadison determinant. For non- invertible X, the multiplicativity 
property sometimes fails. However, it holds if a certain condition on images and 
domains of the multiplicands is fulfilled: 

Proposition 1 Let V be the closure of the range of the operator X. If Y is 
an infective mapping on V and is the zero operator on V , then det (YX) — 
det (Y) det (X) . 

The claim of this proposition is a direct consequence of Theorem 3.14 and 
Lemma 3.15(7) in @]. 

3 Definition of Lyapunov exponents for free op- 
erators 

A pair (A, E) is a tracial W*- non- commutative probability space if A is a finite 
von Neumann algebra with a normal faithful tracial state E, and E (I) = X. 
The trace E will be called the expectation by analogy with classical probability 
theory. 

Let Ai....,A n be sub- algebras of algebra A, and let a* be elements of these 
sub-algebras such that at £ Am) . 

Definition 2 The sub-algebras Ai t ...,A n (and their elements) are called free 
or freely independent if E {a\...a m ) = whenever the following two conditions 
hold: 

(a) E (af) — for every i, and 

(b) k(i) ^ k (i + 1) for every i < m. 



Mx (S) = E(1 S (X)), 
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The random variables are called free or freely independent if the algebras 
that they generate are free. (See 0, or 
of free probability theory.) 

Let {Xi}°Z 1 be a sequence of free identically-distributed operators. Let 
n„ = X n . . . X\, and let Pt be a projection which is free of all Xi and has the 
dimension t, i.e., E (P t ) = t. 

Definition 3 The integrated Lyapunov exponent corresponding to the sequence 
Xi is a real-valued function of t £ [0, 1] which is defined as follows: 

F(t) = lim -logdet(II„Pt). 

n — >oo n 

provided that the limit exists. 

Remark: In the case of random matrices, ±1„ is the product of independent 
identically-distributed random matrices. In this case, it turns out that the 
function defined analogously to F (t) equals the sum of the tN largest Lyapunov 
exponents divided by TV, where TV is the dimension of the matrices and t belongs 
to the set {0/iV, 1/N,..., N/N} . 

Our first task is to prove the existence of the limit in the previous definition. 

Theorem 1 Suppose that Xi are free identically- distributed operators in a tra- 
cial W* -probability space A with trace E. Let u =: dimker (Xi) . Then 

F(t) = i 3 limAio Slog +A {PtXiXiPt) , ift<l-u, 
U \ \ lim A10 E log +A {P U XIX X P U ) , ift > 1 - u. 

Before proving this theorem, let us make some remarks. First, this theorem 
shows that the integrated Lyapunov exponent of the sequence {Xi} exists and 
depends only on the spectral distribution of X*X { . 

Next, suppose that we know that F (t) is differentiable almost everywhere. 
Then we can define the marginal Lyapunov exponent as / (t) — F' (t) . We can 
also define the distribution function of Lyapunov exponents by the formula: 
T(x) = fj, {t € [0, 1] : / (t) < x} , where [i is the usual Borel-Lebesgue measure. 
Intuitively, this function gives a measure of the set of the Lyapunov exponents 
which are less than a given threshold, x. In the finite-dimensional case it is 
simply the empirical distribution function of the Lyapunov exponents, i.e., the 
fraction of Lyapunov exponents that fall below the threshold x. 

Proof of Theorem [TJ The proof is through a sequence of lemmas. We 
will consider first the case of injective operators Xi and then will show how to 
generalize the argument to the case of arbitrary Xi. 

Let Pa denote the projection on the closure of the range of operator A. 
In the following lemmas we always assume that operators belong to a tracial 
W^-probability space A with trace E. 
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Lemma 1 Suppose that operator A is injective, and that P t is a projection of 
dimension t. Then projection Pap ± is equivalent to Pt- In particular, E (Pap 4 ) = 
t. 

Proof: Recall that polar decomposition is possible in A. (See Proposition 
II. 3. 14 on p. 77 in [16[ for details.) Therefore, we can write AP t = WB, where 
W is a partial isometry and B is positive, and where both W and B belong 
to A. By definition, the range of W is [Range (AP t )] , and the domain of W 
is [x : Bx = 0} 1 - = [x : AP t x = Q] X = [x : P t x = 0]^ = [Range (P t )] . Therefore, 
PAP t is equivalent to Pt, with the equivalence given by the partial isometry W. 
In particular, dim (Pap ( ) = dim (P t ) , i.e., E {Pap± ) = t- QED. 

Lemma 2 If A, A*, and P t are free from an operator subalgebra B, then Pap ± 
is free from B. 

Proof: PAP t belongs to the W*-algebra generated by /, A, A* , and Pt- By 
assumption, this algebra is free from B. Hence, Pap ± is also free from B. QED. 

Let us use the notation Q k = Px k ...x 1 p t for k > 1 and Qo = Pt- Then 
by Lemma [U Qk is free from Xk+i- Besides, if all Xj are injective, then their 
product is injective and, therefore, by Lemma [T] Qk is equivalent to Pt. 

Lemma 3 If all Xi are injective, then 

nn 
det {XtQi-x). 
%—\ 

Proof: Note that fI n P t = X n Q n _\X n _\ . . .QiXiQq. We will proceed by 
induction. We need only to prove that 

det (X k+1 Q k X k . . . Q1X1Q0) = det {X k+1 Q k ) det (X k . . . QiXxQo) . (6) 

Let V k be the closure of the range of X k . . . QiXiQq. Since X k+ i is injective 
and Q k is the projector on V k , therefore X k+ iQ k is injective on V k and equal to 
zero on V k . Consequently, we can apply Proposition [1] and obtain ©. QED. 

Now we are ready to prove Theorem [T] for the case of injective Xi. Using 
Lemma [3l we write 

i 1 x — ^ n 

nr 1 log det (n„P t ) = - V log det (XiQi-i) . 
n z — 'i=l 

Note that Xi are identically distributed by assumption, Qi have the same di- 
mension by LemmaQ] and Xi and Qi-% are free by LemmaO This implies that 
lim^o £'log +A {Qi-\X* XiQi-i) does not depend on i, and hence, det 
does not depend on i. Hence, using i = 1 we can write: 

n- 1 log det (n„P t ) = log det (XiP t ) . 

This finishes the proof for the case of injective Xi. For the case of non- 
injective Xi, i.e., for the case when dimker(Xi) > 0, we need the following 
lemma. 
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Lemma 4 Suppose that P t is a projection operator free of A and such that 
E(P t ) =t. Then dim ker(AP t ) = max{l - t, dimker {A)} . 



Proof: Let V = (Ker A) and let Py be the projection on V. Then E (Py) = 
1 — dim Ker A. Note that Ax = <^=> Pyx — 0. Consequently, AP t x = <^=> 
PyP t x = 0. Therefore, we have: 

dim {x : AP t x = 0} = dim {x : PyP t x = 0} 

= dim {x : P t PyP t x = 0} . 

Since Pt and Py are free, an explicit calculation of the distribution of PtPyPt 
shows that 

dim {x : P t PyP t x = 0} = max {1 - t, 1 - dim V} . 

QED. 

Consider first the case when < dim ker Xi < 1 — t. This case is very similar 
to the case of injective Xi. Using Lemma|4]we conclude that dim Ker (X\Pt) = 
1 — t, and therefore that E (Px 1 p t ) = t. If, as before, we denote PxiP t as Qi, 
then the projection Q\ is free from X2, and E (Qi) = t. 

Similarly, we obtain that E (Px 2 Qt) = t. Proceeding inductively, we define 
Qk = Px h Q k _ 1 and conclude that Qk is free from Xk+i and that E (Qk) = t. 

Next, we write X k ...XiP t = X fc Q fe _ 1 X fc _ 1 Q fc _ 2 ---^i ( 3o; where Q denotes 
P t , and note that X^Qk-i is injective on the range of Qk-i- Indeed, if it were 
not injective, then we would have dim (Ker (Xf.) H Range (Qk-x)) > 0- But this 
would imply that dim (Ker XkQk-i) > dim (KerQfe-i) = X — t, which con- 
tradicts the fact that dim (Ken XkQk-i) = X — t. Therefore, Proposition Q] is 
applicable and 

6et(X k ...X 1 P t ) = det(X fc Q fc _ 1 )...det(XiQo) 
= [dot {X 1 P t )] k . 

Now let us turn to the case when dim ker Xj = u > X—t. Then dim Ker (XiP t ) 
1 and therefore E (Px 1 p t ) — u. Proceeding as before, we conclude that 
E(Qk) =MforallA:> 1, and we can write Xk—XxP t = X k Qk-iXk-iQk-2---XiQ , 
where we have denoted Pt as Q$. Then we get the following formula: 

det{X k ...X x P t ) = det(X fe Q fc _ 1 )...det(X 2 Q 1 )det(X 1 Q ) 
= [det (XiP„)] &_1 det (XiPt) . 

Therefore, 

lim n^ 1 log det (II n P t ) = log det (X 1 P U ) . 

n^oo 

QED. 
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4 Example 



Let us compute the Lyapunov exponents for a random variable X that has the 
product X*X distributed according to the Marchenko-Pastur distribution. Re- 
call that the continuous part of the Marchenko-Pastur probability distribution 



with parameter A > is supported on the interval 
and has the following density there: 



i-Va) , (i + Va)' 



P\ (x) = — 



2nx 

For A 6 (0, 1) , this distribution also has an atom at with the probability mass 
(1 — A) assigned to it. The Marchenko-Pastur distribution is sometimes called 
the free Poisson distribution since it arises as a limit of free additive convolutions 
of the Bernoulli distribution, and a similar limit in the classical case equals the 
Poisson distribution. It can also be thought of as a scaled limit of the eigenvalue 
distribution of Wishart-distributed random matrices (see for a discussion). 

Proposition 2 Suppose that X is a non- commutative random variable in a fra- 
ctal W* -probability space (A, E) , such that X* X is distributed according to the 
Marchenko-Pastur distribution with parameter A. If A > 1, then the distribution 
of Lyapunov exponents of X is 

f 0, tfa;<(l/2)log(A-l) 
F(x) = l e 2x + l-X, ifxe [±log(A- l),±logA) . 
I 1 i/z>±logA. 

If X < 1, then the distribution of Lyapunov exponents of X is 

( e 2 *, ifx< (1/2) log (A) 
T{x) = l A, ifx£ [ilog(A),0) 
I 1 ifx>0. 

Remark: If A = 1, then the distribution is the exponential law discovered 
by C. M. Newman as a scaling limit of Lyapunov exponents of large random 
matrices. (See 10 1, and This law is often called the "triangle" law 
since it implies that the exponentials of Lyapunov exponents converge to the 
law whose density function is in the form of a triangle.) 

Proof of Proposition [2) It is easy to calculate that the continuous part of 



the distribution olP t XP t is supported on the interval 
and has the density function 

y/4\t-[x-(t + A)] 2 
Pt,x \x) = 



(Vi-Vx) , (Vi + \/A)' 



2irx 
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This distribution also has an atom at x = with the probability mass max {1 — A, 1 
See for example, results in [12 1. 

Next, we write the expression for the integrated Lyapunov exponent. If 
A > 1, or A < 1 but A > t, then 

1 



Fx(t) 



-\imE\og +s (P t X*XP t ) 

2 e|0 



(Vt+\/X) 2 



logx 



4At -[ x -(t + X)Y 



2irx 



dx. 



(7) 



If A < 1 and A < t, then 



1 



F x (t) = ilim^log +e (P 1 _ A X*XP 1 _ A ) 

Z elO 



^4A(l-A)-[:r-l] 2 
log a;- 1 ax. 



2irx 



(8) 



Differentiating ([7]) with respect to t, we obtain an expression for the marginal 
Lyapunov exponent: 



h (t) = 



l 

-in 



(v-t+vT) 2 



log a; 



x - t + A 



(Vt-v'X) 2 



i eta;. 



4Air - [x - t + A] z 



(9) 



/ (2VAt) - 1 and then 9 



Using substitutions u — 
arccosu, this integral can be computed as 

/ A (t) = ilog(A-t). 

From this expression, we calculate the distribution of Lyapunov exponents 
for the case when A > 1: 

r 0, if x < (1/2) log (A- 1) 

F{x)=l e 2x + l-\, if 16 [|log(A- 1), I log A) 
1 1 ifx>ilogA. 

A similar analysis shows that for A < 1, the distribution is as follows: 

e 2x , if x< (1/2) log (A) 
T(x) = { A, if.xe [ilog(A),0) 
1 if x > 0. 



QED. 



10 



5 A relation with the ^-transform 



In this section we derive a formula that makes the calculation of Lyapunov 
exponents easier and relates them to the S'-transform of the operator Xi . Recall 
that the ■(/'-function of a bounded non-negative operator A is defined as ip A (z) = 
J2T=i E i Ak ) zk - Then the ^-transform is S A (z) = (l + z _1 ) V4~ X) (z) , where 
tp A ^ (z) is the functional inverse of ipA ( z ) m a neighborhood of 0. 

Theorem 2 Let Xi be identically distributed free bounded operators in a tra- 
cial W* -probability space A with trace E. Let Y — X^X\ and suppose that 
/i Y ({0}) = 1 — u € [0, 1) , where fi Y denotes the spectral probability measure of 
Y. Then the marginal Lyapunov exponent of the sequence {Xi} is given by the 
following formula: 

-|log[5rH)] ift<u, 
ift>u, 



fx (t) = 



where Sy is the S -transform of the variable Y. 

Remark: Note that if Xj* X\ has no atom at zero then the formula is simply 

/*(*) = -! log [M-t)]. 

Proof: If t > u, then fx (i) = by Theorem [T] Assume in the following 
that t < u. Then PtX*XPt has an atom of mass 1 — t at 0. Let \i t denote the 
spectral probability measure of PtX*XPt, with the atom at removed. (So the 
total mass of [x t is t.) We start with the formula: 

logx = \og{c + x) - { 

Jo x + s 

and write: 

/ \ogx fj, t (dx) = lim < t log (c) + G t (— s) ds > , 
Jo c ^°° I Jo J 

where Gt is the Cauchy transform of the measure \i t . 

Next, note that Gt (—s) = —s tp t (— s _1 ) — ts^ 1 and substitute this into 
the previous equation: 



I 



r°° \ r^t( s_1 ) 

' log a; n t (dx) = lim < t log c — t log c + 1 log (e) 4- / — -ds 

c^oc.e^O J e S 



Using substitutions v = — logs, and A — — logs, we can re- write this equation 
as follows: 

f f A \ 

\ogx fx t (dx) = lim < — tA — I ip t (—e v ) dv > . 
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The function ip t (—e v ) monotonically decreases when v changes from — oo to 
oo, and its value changes from to —t. Let s* =: tp t (— e°) = tp t (—1) and let 
£ t (x) denote the functional inverse of tp t (—e v ) . The function £ t (x) is defined 
on the interval (— i, 0) . In this interval it is monotonically decreasing from oo 
to — oo. The only zero of £ t (x) is at x — s* . 

It is easy to see that 

lirn^ S^-tA - jf Vt (-e") dwj = - y* (s) da:, 
- /" ^(-e")di; = - f\ t (x)dx. 

J —oo J s* 



and that 



Therefore, 



/ log ^ Mt (da;) = - / £ t (a;) da;. 
Jo J-t 

It remains to note that £ t (x) = log ^ (x) 

) 

log x /U t (dx) = 



in order to conclude that 
dx. 



log 1} (x) 



The next step is to use Voiculescu's multiplication theorem and write: ip\ ^ (x) 
tp ( Y 1} ( x ) (1 + x) / (t + x) . Then we have the formula: 



f 

Jo 



logx \x t (dx) 



L 



log 



= - / log 



dx 



/: 



log 



1 + x 



i + x 



dx 



1} (x)j dx + (1 - t) log (1 - t) + ilogi. 



The integrated Lyapunov exponent is one half of this expression, and we can 
obtain the marginal Lyapunov exponent by differentiating over t: 



/(*) 



l - (- log [-^ 1} (-*)] + log* - log (1 - t)) 



= -2 log 



1 - I) ("*) 



= -_log[S y (-t)]. 

QED. 
Example 

Let us consider again the case of identically distributed free such that 
X*Xi has the Marchenko-Pastur distribution with the parameter A > 1. In this 
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case Sy (z) = (A + z) 1 . Hence, applying Theorem [51 we immediately obtain a 
formula for the marginal Lyapunov exponent: 



/(t) = -log(A-t). 



Inverting this formula, we obtain the formula for the distribution of Lyapunov 
exponents: 



which is exactly the formula that we obtained earlier by a direct calculation 
from definitions. It is easy to check that a similar agreement holds also for 



In the following corollaries we always assume that operators belong to a 
tracial M / *-probability space A with trace E. 

Corollary 1 Let X and Y be free and such that X* X and Y*Y are bounded 
and have no atom at zero. Let fx, /y, and fxy denote the marginal Lyapunov 
exponents corresponding to variables X, Y and XY, respectively. Then 



Corollary 2 If X is bounded and X*X has no atom at zero, then the marginal 
Lyapunov exponent is (weakly) decreasing in t, i.e. f' x (t) < 0. 

Proof: Because of Theorem [2] we need only to check that S (t) is (weakly) 
decreasing on the interval [—1,0] , and this was proved by Bercovici and Voiculescu 
in Proposition 3.1 on page 225 of [if. QED. 

Corollary 3 If X is bounded and X*X has no atom at zero, then the largest 
Lyapunov exponent equals (l/2)logE(X*X). 




A < 1. 



fxY (t) = fx (t) + fv (t) . 



Proof: By Theorem 



fxY (t) 



--log [Sy-X'XY (-*)] 

-llog[S Y , Y (-t) S x * x (-t)] 
fx (t) + fv (t) . 



QED. 
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Proof: This follows from the previous Corollary and the fact that Sy (0) = 
1/E(Y) . QED. 

Remark: It is interesting to compare this result with the result in Q , which 
shows that the norm of the product of N x N i.i.d. random matrices X\, . . . , X n 
grows exponentially when n increases, and that the asymptotic growth rate ap- 
proaches \ log-E(tr [X\X\)) if N — ► oo and matrices are scaled appropriately. 
The assumption in [2] about the distribution of matrix entries is that the dis- 
tribution of X*Xi is invariant relative to orthogonal rotations of the ambient 
space. Since the growth rate of the norm of the product X\ . . . X n is another 
way to define the largest Lyapunov exponent of the sequence Xi , therefore the 
result in Q is in agreement with Corollary [3l 

The main result of Theorem [2] can also be reformulated as the following 
interesting identity: 

Corollary 4 IfY is bounded, self-adjoint, and positive, and if {P t } is a family 
of projections which are free ofY and such that E (P t ) = t, then 



log Sy(-t) 



d 

~dt 

= -4 

dt 



lim £;iog +A (P t YP t ) 
logdet (VYP t 



Conversely, we can express the determinant in terms of the S'-transform: 
Corollary 5 If X is bounded and invertible, then 



det(X) =expj-i^ \ogS x * 



x 



(-*) dt 



6 Conclusion 



In conclusion we want to indicate how our results are related to results by New- 
man in Q. Newman considers N-by-N random matrices Xi with rotationally- 
invariant distribution. of entries. Suppose that the empirical distribution of 
eigenvalues of yjX*Xi converges to a probability measure K (dx) as N — > oo. 
Let H [x) be the limit cumulative distribution function for e^ k , where ft are 
Lyapunov exponents of Xi. Newman shows that this limit exists and satisfies 
the following integral equation: 

°° 2 

K{di) = \. 



o 



H{x) x 2 + (1-H(x))t 



2 - 



We are going to transform this expression in a form which is more compa- 
rable with our results. Let fi (dt) be the limit of the empirical distribution of 
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eigenvalues of X*Xi as N — > oo, and let s — t 2 . Then, we can re- write Newman's 
formula as follows: 



/i (ds) = 1. 



if (x) x 2 + (1 - if (x))s 



This expression can be transformed to the following equality: 

where i/v (2) is the i/'-fimction f a non-commutative random variable Y with 
the spectral probability distribution /i. Hence 

^ = S Y [H (x) - 1] . 

If we denote 1 — if (x) as u, then 

X = e - 2 lo 8[-Sv (-«)]_ 

This implies that as N grows, the ordered Lyapunov exponents ff. converge 
to — j log [Sy (— it)] , provided that fe/iV — > u. This is in agreement with our 
formula in Theorem O 
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